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Abstract 

Background: FeFe-hydrogenases are the most active class of H2-producing enzymes known in nature and may 
have important applications in clean H2 energy production. Many potential uses are currently complicated by a 
crucial weakness: the active sites of all known FeFe-hydrogenases are irreversibly inactivated by O2. 

Results: We have developed a synthetic metabolic pathway in £ coli that links FeFe-hydrogenase activity to the 
production of the essential amino acid cysteine. Our design includes a complementary host strain whose 
endogenous redox pool is insulated from the synthetic metabolic pathway. Host viability on a selective medium 
requires hydrogenase expression, and moderate O2 levels eliminate growth. This pathway forms the basis for a 
genetic selection for O2 tolerance. Genetically selected hydrogenases did not show improved stability in O2 and in 
many cases had lost H2 production activity. The isolated mutations cluster significantly on charged surface residues, 
suggesting the evolution of binding surfaces that may accelerate hydrogenase electron transfer. 

Conclusions: Rational design can optimize a fully heterologous three-component pathway to provide an essential 
metabolic flux while remaining insulated from the endogenous redox pool. We have developed a number of 
convenient In vivo assays to aid in the engineering of synthetic H2 metabolism. Our results also indicate a H2- 
independent redox activity in three different FeFe-hydrogenases, with implications for the future directed evolution 
of H2-activating catalysts. 



Background 

Nature provides a catalyst for H2 production: hydroge- 
nase metalloenzymes[l]. Hydrogenase active sites are 
composed only of widely available metals, yet can 
reduce protons to H2 more rapidly than platinum cata- 
lysts[2]. They have therefore been proposed as a substi- 
tute for precious metal electrodes for the low-cost, high- 
volume interconversion of electricity and H2[2,3]. Enzy- 
matic catalysts are particularly well-suited for synthetic 
integration with specifically biological energy sources. 
For example, coupling of hydrogenase enzymes with the 
photosynthetic machinery enables the direct production 
of H2 from sunlight, as has been demonstrated both in 
vitro[4:] and in vivo[5,6]. 
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Several classes of hydrogenases are found in nature[l], 
distinguished primarily by the metal content of their 
active site. These enzymes vary widely with respect to 
their activity level, maturation requirements, redox part- 
ners and O2 tolerance, all of which affect their potential 
biotechnological utility. 

Here we employ iron-iron (FeFe) hydrogenases of the 
monomeric, cytosolic class found widely in nature and 
common among the Clostridia[7]. These are the bio- 
chemically best-characterized hydrogenases, featuring a 
well-defined set of maturation factors that are sufficient 
for their heterologous expression in engineered organ- 
isms [8,9]. Members of this class exhibit among the high- 
est activities of known hydrogenases, with catalytic rates 
exceeding the 1000 nmol H2 mg' min'^ estimated to be 
required for efficient in vivo photosynthetic H2 produc- 
tion[10]. However, biotechnological applications of 
FeFe-hydrogenases are impeded by the rapid and irre- 
versible inactivation of the active site by 02[11]. 
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In contrast, several naturally occurring hydrogenases 
of the nickel-iron (NiFe) class[12-15] and the nickel- 
iron-selenium (NiFeSe) class[16] are known to be O2- 
tolerant. However, the H2 production rates of NiFe- 
hydrogenases are 10 to 1000 times slower than FeFe- 
hydrogenases[l,17,18]. The maturation requirements for 
NiFe-hydrogenases are also relatively complex, with as 
many as 13 genes required for heterologous[l,17,18] 
reconstitution in some cases [19], limiting their facility of 
heterologous expression. 

There is currently no known hydrogenase of any class 
that produces H2 at rates comparable to the FeFe- 
hydrogenases and is tolerant of O2. However, both the 
activities and O2 tolerance of naturally occurring FeFe- 
hydrogenases span wide ranges and are uncorrelated 
[20]. O2 sensitivity may even be adaptive in hydrogenase 
evolution, acting as a form of regulation[21,22]. If the 
natural selective pressure for O2 tolerance is weak, it 
should be possible to engineer a hydrogenase with both 
high activity and improved O2 tolerance. 

Recent results suggest structural modifications that 
may alter hydrogenase O2 tolerance. Widening a gas 
channel, through which H2 is believed to approach the 
active site, resulted in the loss of O2 tolerance in a Ral- 
stonia eutropha NiFe-hydrogenase[23]. Similarly, the 
mutation of gas channel residues in a Desulfovibrio fruc- 
tosovorans NiFe-hydrogenase increases O2 tolerance yet 
only moderately reduces catalytic activity[24] . Although 
understanding of the mechanism of 02-inactivation of 
hydrogenases has greatly improved [25,26] it has not yet 
been possible to rationally design mutations to improve 
O2 tolerance without impairing activity. 

Directed evolution permits the optimization of func- 
tion without detailed knowledge of enzyme mechanism 
[27,28]. This technique entails either screening or 
genetically selecting for desirable mutants from suitably 
large mutation libraries. In vitro screens for hydrogenase 
activity are capable of assaying 10 - 10 hydrogenase 
mutants per day [29]. Yet double or triple combinatorial 
point mutation libraries for a gene the size of the hydro- 
genase approach 10^ - 10^'' unique sequences, and can 
not be comprehensively screened with existing methods. 
Genetic selection is the most powerful available tool for 
the rapid and economical parsing of very large mutation 
libraries. This technique requires a connection between 
a desired enzyme property and the evolutionary fitness 
of a strain expressing the enzyme. 

We have engineered an artificial pathway in E. coli 
connecting hydrogenase function to the production of 
sulfide, an essential precursor of the amino acids 
cysteine and methionine. We address two design objec- 
tives that maximize the utility of this system for directed 
evolution. First, we minimize pathway activity indepen- 
dent of the hydrogenase, insulating redox interactions 



with native metabolism and reducing the potential for 
selection false positives. Second, we maximize hydroge- 
nase-dependent activity of the pathway, creating a 
robust link to host fitness by ensuring binding compat- 
ibility between synthetic pathway components. Because 
genetic selection can be performed only in vivo, we have 
optimized the biochemistry of our pathway in vivo. This 
pathway forms the basis of a genetic selection for O2- 
tolerant hydrogenases. 

Methods 

Cloning and gene synthesis 

All cloning was performed in E. coli DH5a using stan- 
dard BioBrick assembly techniques [30]. Final constructs 
were assembled in commercial Duet vectors (Novagen) 
with multiple cloning sites modified to accept BioBrick 
parts. The plasmids used in these experiments are listed 
in table 1. Complete vector sequences are provided in 
additional file 1. 

The hydrogenase genes hydAl from Clostridium acet- 
obutylicum ATCC 824 and hydA from Clostridium sac- 
charobutylicum P262 were cloned from plasmids 
received from Matthew Posewitz (National Renewable 
Energy Laboratory, Golden, CO, USA). Dr. Posewitz 
also provided plasmids bearing the C. acetobutylicum 
maturation factors hydE, hydF, and hydG['il]. The 
hydrogenase maturation factors hydEF and hydG from 
Chlamydomonas reinhardtii were commercially synthe- 
sized by Codon Devices. The par gene encoding pyru- 
vate-ferredoxin oxidoreductase from Desulfovibrio 
africanus was cloned from plasmid pLPl[32] provided 
by Laetitia Pieulle (Centre National de la Recherche 
Scientifique, Marseille, France). Clostridium acetobutyli- 
cum ferredoxin CAC0303[33] was cloned from genomic 
DNA (ATCC 824). Spinacia oleracea ferredoxin (Gen- 
Bank AAA34028), Zea mays ferredoxin (GenBank 
ACA34367) and Zea mays sulfite reductase (GenBank 
BAA23641) were cloned from total RNA. The ferre- 
doxin-NADPH reductase (FNR) gene from Zea mays 
(GenBank AAB40034) was synthesized by Codon 
Devices. Chloroplast transit peptides were omitted from 
all plant-derived constructs. The codon usage of syn- 
thetic genes was optimized by the manufacturer for het- 
erologous expression. 

Gene expression 

All synthetic pathways were expressed in strains derived 
from E. coli BL21(DE3), obtained from Agilent Technol- 
ogies. Plasmids were maintained through selective anti- 
biotics: ampicillin, 50 i^g/mL; spectinomycin, 25 |ig/mL; 
kanamycin, 25 |ig/mL; chloramphenicol, 12.5 |ig/mL. 
Protein expression from the Duet vectors was T7-pro- 
moter driven and induced with IPTG at a concentration 
of 1 mM. 
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Table 1 Plasmids used in this study 



Name 


Constructs 


Backbone 


Resistance 


Source 


Hydrogenase activity in wVo* 


„ 1 — r „ 1 

pt 1 .mpl 


caHydt coHydA 


pETDuet-1 


Ampicillin 


Matthew Posewitz[31] 


pCDF.mp2 


caHydF caHydC 


pCDFDuet-1 


Spectinomycin 


Matthew Posewitz[31] 


pACYCewB 


doPFOR 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYC.ew4 


doPFOR soFD 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewS 


doPFOR zmFD 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYC.ew6 


daPFOR ctFD 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCew? 


doPFOR cam 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewS 


5oFD 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYC.ew9 


zmFD 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewl 0 


cfFD 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewl 1 


cam 


pACYCDuet-1 


Chloramphenicol 


This work 


FNR-supported growth^ 


pCDF.ewl2 


zmFNR 


pCDFDuet-1 


Spectinomycin 


This work 


pACYCewl 3 


5oFD zmSIR 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewl 4 


zmFD zmSIR 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewl 5 


cfFD zmSIR 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewl 6 


cam zmSIR 


pACYCDuet-1 


Chloramphenicol 


This work 


pACYCewl? 


zmSIR 


pACYCDuet-1 


Chloramphenicol 


This work 


Hydrogenase-supported growth and selection*" 


pACYCewl 8 


crHydEF crHydG 


pACYCDuet-1 


Chloramphenicol 


This work 


pET.ew19 


5oFD zmSIR 


pETDuet-1 


Ampicillin 


This work 


pCDF.ew20 


cfHydA 


pCDFDuet-1 


Spectinomycin 


This work 


pCDF.ew21 


caHydA 


pCDFDuet-1 


Spectinomycin 


This work 


pCDF.ew22 


csHydA 


pCDFDuet-1 


Spectinomycin 


This work 



Duet vector backbones (Novagen) were used for all protein expression. Complete vector sequences are provided in additional file 1. A) Plasmids mpl-ewl 1 were 
used to generate PFOR-driven hydrogenase activity data for Figure 4. B) Plasmids ewl2-ew17 were used to assess knockout strain insulation for Figure 2. They 
were also used to measure ferredoxin performance in Figure 4. C} Plasmids ew18-ew22 facilitated the hydrogenase O2 tolerance measurements in Figure 3. They 
were also used in the hydrogenase-supported growth curves in Figure 5 and for the genetic selection. 



Gene deletions 

Sequential gene deletions were constructed by PI phage 
transduction from the Keio collection [34]. Serial dele- 
tions were enabled by removing the transduced kanamy- 
cin resistance marker though homologous recombination 
at flanking FRT sites. Transient expression of the Flp 
recombinase was facilitated by the 705-Flp plasmid, 
which exhibits temperature-sensitive recombinase 
expression and a temperature-sensitive replication origin 
[35]. Kanamycin marker integration and subsequent 
removal was confirmed for all loci by PCR. 

Selective and Induction media 

Selective media was a standard M9 formulation, supple- 
mented with additional glucose, sulfate, ferric iron and a 
rich mix of supplements less cysteine and methionine. 
Induction media for hydrogenase expression was LB 
with added glucose, ferric iron, phosphate buffer and 
Baker's antifoam reagent. Exact media recipes are pro- 
vided in additional file 2. 



Anaerobic technique and custom atmospheres 

Anaerobic liquid culture was performed in 40 mL 
serum vials sparged with N2 and sealed with SubaSeal® 
rubber septa (Sigma-Aldrich). To maintain anaerobio- 
sis during handling, samples were drawn and reagents 
added by piercing the septa with non-coring syringe 
needles. 

Agar plates were incubated under defined gas mix- 
tures within sealed Vacu-Quick jars (Almore Interna- 
tional). The ambient atmosphere was removed by 
several cycles of evacuation and replacement with pure 
N2 before supplying a custom atmosphere. Aluminosili- 
cate desiccant packets were added to prevent moisture 
accumulation within the jars. 

Growth assays 

Cells were grown to saturation in induction media and 
washed 3x with phosphate-buffered saline (PBS). Cells 
were resuspended in selective media at an initial ODgoo 
of 0.01. Final ODs were measured after 18 hours of 
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growth at 37°C. Anaerobic conditions, when appropriate, 
were introduced as described above. 

In situ hydrogenase activity assays 

Cells were grown to saturation in anaerobic induction 
media. Samples were drawn to determine cell density by 
OD. Fresh serum vials containing 25 mL of induction 
LB were anaerobically inoculated with 10 cells (~ 5 
mL). Following incubation at 37°C for 2 hours, H2 pro- 
duction was stopped by the addition 2 mL methanol. 
Accumulated headspace H2 was measured by gas chro- 
matography (Shimadzu GC-14A). 

In vitro hydrogenase activity assays 

Hydrogenase activities were measured with a biochem- 
ical methyl viologen assay adapted from King et al.[9]. 
Hydrogenase-expressing E. coli were grown to satura- 
tion in 20 mL of induction media under anaerobic 
conditions. Samples were drawn to determine cell den- 
sity by OD. Cells were lysed with 1 mL of lysis buffer 
consisting of 20 mL B-PER II protein extraction 
reagent (Thermo Scientific), 500 \iL Baker's antifoam 
B, 100 units DNAse I and 50 mg dithionite. Antifoam 
and DNAse I were added to prevent foaming of the 
lysate. Dithionite served to scavenge dissolved O2 in 
the buffer. Lysis continued for 15 minutes under con- 
tinuous N2 sparging. Lysed cultures were sealed and 
injected with 1 mL of methyl viologen assay buffer 
consisting of 20 mL 1 M Tris at pH 8, 300 mg methyl 
viologen and 3 g dithionite. Following 2 hours of incu- 
bation at 37°C, H2 production was stopped with 2 mL 
of methanol. Accumulated headspace H2 was measured 
by gas chromatography. 

Sulfide production assays 

Sulfide production was measured spectrophotometrically 
by the methylene blue method [36]. Bacteria were grown 
to saturation in selective media with IPTG and appro- 
priate antibiotics, supplemented with 100 mg/L cysteine 
to allow growth of all strains. Cultures were diluted 1:50 
to 25 ml in 40 ml sealed serum vials. Vials were flushed 
with either pure N2 or 40% H2 in N2 and grown for 6 
hours at 37°C. Cuvettes were prepared with 2.5 mL 
assay solution containing N,N-Dimethyl-l,4-phenylene- 
diammonium dichloride (200 |iM), FeCl (600 \xM) and 
HCl (60 mM). 25 \iL of filtered media was added to 
assay cuvettes and allowed to react for 1 hour. Absor- 
bance was measured at 660 nm and compared to a stan- 
dard curve to calculate sulfide concentrations. Statistics 
were inferred from three biological replicates. 

Hydrogenase O2 tolerance assays 

Hydrogenase half-lives in O2 were measured using an 
adaptation of the in vitro activity assay above. Following 



anaerobic lysis, a customized sparging apparatus was 
used to bubble O2 through the cultures at a constant 
partial pressure of 0.1 atm for defined time points 
between 5 and 25 minutes. Cultures were then flushed 
for 5 minutes with pure nitrogen before being sealed 
and assayed with methyl viologen assay buffer, as above. 
Measurements at each time point were taken for 3 bio- 
logical replicates. 

Colony growth assays under custom atmospheres 

Cells were grown to saturation in induction media. Resi- 
dual nutrients were removed by 3x washing with PBS. 
Cell densities were measured with a hemocytometer and 
diluted to a final concentration of 1 cell/^L. Fifty micro- 
liters (~ 50 cells) were dispensed onto selective plates 
and transferred to Vacu-Quick jars. The jars were filled 
with 15% H2, O2 varying from 0-12.5%, and a balance of 
N2 to a total internal pressure of 1 atm. Incubation was 
carried out at 37°C for 72 hours. Plates were photo- 
graphed in an inverted camera stage (Figure 5E). Colo- 
nies were identified and sized with an image analysis 
script implemented in MATLAB (MathWorks). Each 
data point represents size data collected from roughly 
50 individual colonies. 

Hydrogenase mutagenesis and selection 

Mutation libraries were created through error-prone 
PCR amplification of the hydrogenase genes with the 
GeneMorph II Random Mutagenesis Kit (Agilent Tech- 
nologies). Exact primer sequences are indicated in addi- 
tional file 1. The PCR products were digested with NotI 
and Spel (for caHydA) or Ncol and NotI (for csHydA), 
purified, and ligated into a modified pCDFDuet-1 vector 
(Novagen). A target frequency of 4.5 mutations/kb was 
sought according to the manufacturer's instructions and 
confirmed by sequencing several clones plated 
nonselectively. 

Host cells were transformed with the mutant hydroge- 
nase libraries by electroporation and recovered in SOC 
media for 1 hour at 37°C. Following recovery, the cells 
were washed three times with Phosphate Buffered Saline 
(PBS) and plated on selective media. Approximately 10 
unique transformants were applied to each plate, as 
assessed by serial dilution of the recovered cells. Selec- 
tion plates were then transferred to Vacu-Quick jars 
prepared as above with a selective atmosphere contain- 
ing 10% O2 and 15% H2. Incubation proceeded at 37°C 
for 4 days. 

Selection restreaks were performed using Petri dishes 
with internal divisions to block the diffusion of metabo- 
lites from neighboring streaks [37]. Soft plastic inoculat- 
ing loops were used to prevent marring of the agar 
surface and maintain a uniformly oxic environment 
[38,39]. 
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Structure and surface modeling 

Homology models of hydrogenase mutants were con- 
structed using SWISS-Model homology modeling ser- 
vice[40] using the Clostridium pasteurianum 
hydrogenase I structure[41] as a template (PDB ID: 
3C8Y). 

The protein surface was inferred with the NIH MBI 
Diffusion Accessibility Calculator[42]. The diffusion 
accessibility was averaged in a moving window of 5 resi- 
dues, then thresholded to produce defined surface- 
accessible regions. 

Simulated mutation spectra 

The set of point mutations isolated by our genetic selec- 
tion was compared to the expected distribution under a 
null model of no selective pressure. The null distribu- 
tion was generated as a parallel Monte Carlo simulation 
implemented in PYTHON. Our algorithm produced 
5-10^ randomly mutated sequences for each caHydA 
and csHydA in batches of 8 and 12, the number of var- 
iants of each enzyme isolated by the genetic selection. 
We used the point mutation spectrum provided by Agi- 
lent for the GeneMorph 11 kit and a total frequency of 
4.5 mutations/kb. Nonsilent mutations were grouped by 
charge of the original and mutant amino acid (i.e. a 
positive to neutral substitution). Statistical 95% confi- 
dence intervals for the null distribution were derived 
from the size of each charge change group across simu- 
lated replicates. 

Results and Discussion 

Design for hydrogenase-dependent growth 

Figure 1 details the synthetic pathway designed for our 
selection and expressed in E. coli. An exogenous FeFe- 
hydrogenase consumes H2 and reduces ferredoxin. Ferre- 
doxin donates electrons to sulfite reductase for the reduc- 
tion of sulfite to sulfide. Sulfide serves in the host as an 
indispensible sulfur source for cysteine biosynthesis. 



6 3 H2O 

H2SO3 V, ^ i^^g 

* *■ — 6 Hp 



12 H 



cysteine, 
methionine, 
other essential 
metabolites 



(HydA) 



Figure 1 A synthetic pathway linking hydrogenase activity to 
host viability. A hydrogenase enzyme (HydA) consumes H2 and 
reduces ferredoxin (FD), a redox carrier protein. Ferredoxin then 
donates electrons to a plant-type sulfite reductase (SIR), which 
reduces sulfite to sulfide. With the native sulfite reductase deleted, 
this pathway becomes the only source of reduced sulfur, an 
essential metabolic precursor 



Ferredoxin, an electron-carrying iron-sulfur (Fe-S) 
protein, is the native redox partner of the best-charac- 
terized FeFe-hydrogenases[43]. A ferredoxin also 
receives electrons from photosystem I in plants, suggest- 
ing that it could be adapted to mediate light-driven H2 
production. We chose ferredoxin as an intermediate 
with this future application in mind. A ferredoxin 
homolog, fdx, is found in E. coli, where it plays an 
essential role as a scaffold site for iron-sulfur cluster 
assembly [44]. E. coli does not appear to use ferredoxin 
as an electron carrier in the metabolic network, instead 
relying on NAD(P)H. We therefore hypothesized that 
ferredoxin chemistry would be insulated from native 
metabolism. 

The KEGG pathway database [45] identifies three 
enzymes that produce essential metabolites using ferre- 
doxin as an electron source. Glutamate synthase, nitrite 
reductase and sulfite reductase activities are all essential 
for the growth of E. coli in minimal medium. The native 
bacterial enzymes draw electrons from NADPH to pro- 
duce glutamate, ammonia and sulfide, respectively. The 
analogous enzymes in plants yield the same products 
while drawing electrons from ferredoxin, a common 
redox carrier in those species. While all three products 
are indispensable for E. coli viability, sulfide is consumed 
in the smallest molar quantity [46]. We therefore chose 
to employ sulfite reductase, reasoning that the small 
metabolic requirement would afford more tolerance for 
suboptimal performance of the heterologous pathway. 

Sulfite reductase is not essential on cysteine-contain- 
ing rich media, but it becomes essential on selective 
media containing only oxidized sulfur sources such as 
sulfate or sulfite. In the absence of the native E. coli 
NADPH-dependent sulfite reductase, cysl, the synthetic 
pathway is the only metabolic source of reduced sulfur. 
If the components of this pathway are insulated from 
any endogenous electron sources, then hydrogenase 
activity will also be essential for growth. Increasing O2 
concentrations, by inactivating the hydrogenase, will 
eventually inhibit the ability of host cells to grow on sul- 
fite. This synthetic pathway therefore enables a genetic 
selection for hydrogenase mutants with an ability to 
support growth in high O2. 

Synthetic ferredoxin-dependent sulfite reduction 

We first sought to establish that the native E. coli BL21 
(DE3) sulfite reductase could be replaced with a ferre- 
doxin-dependent pathway (Figure 2A). Deletion of the 
cysl sulfite reductase did not impair growth on rich 
media (LB), but eliminated growth on selective media 
with sulfate as the sole source of metabolizable sulfur. 
Expression of corn {Zea mays) sulfite reductase (zwSIR) 
alone, or together with spinach ferredoxin (50FD), failed 
to rescue growth. This indicated that the E. coli host 
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B 

Strain 
EW1 
EW2 
EW3 
EW4 
EW5 
EW6 
EW7 
EW11 

c 



Deletions 

cysl 

cyst fpr 

cysl fpr ydbK 

cysl fpr ydbK her 

cysl fpr ydbK her yeaX 

cysl fpr ydbK her yeaX hcaD 

cysl fpr ydbK her yeaX hcaD frdB 

cysl fpr ydbK her yeaX hcaD frdB hycE hyaS hybC hyfG 

■ with FNR 




WT EW1 EW2 EW3 EW4 EW5 EW6 EW7 EW1 1 



Figure 2 Genetic insulation of the synthetic pathway. A) £ co!i 
BL21(DE3) cysl cells were transformed with plasmids expressing 
zmFNR, 5oFD and zmSIR. The cysl deletion conveys a requirement 
for reduced sulfur, which the heterologous pathway supplies. Cells 
were grown for 24 hours on selective media with or without 
atmospheric O2. All three factors were required to rescue growth 
under aerobic conditions (left plate). Expression of zmFNR was not 
required under anaerobic conditions (right plate), indicating that 
50FD was receiving electrons from another source. B) Genetic 
deletions were targeted to eliminate potential endogenous redox 
partners for ferredoxin, therefore linking sulfide production 
specifically to a synthetic electron source. Also deleted were the 
catalytic subunits of each native hydrogenase, ensuring that only 
exogenous H2 was present in our system. C) Each deletion strain 
was transformed with soFD, zmSIR and either zmFNR or an empty 
plasmid. Growth was measured after 18 hours at 37°C under strict 
anaerobic conditions, as described in the methods. Sequentia 
deletions reduced the nonspecific anaerobic background growth, 
with the largest effect produced by the deletion of ydbK. The final 
deletion strain, EW1 1, showed no growth defect in rich media and 
was used in all later experiments. 



provides no interacting source of reduced ferredoxin or 
ferredoxin-reductase activity under these conditions. 

We then provided a heterologous source of reduced 
ferredoxin in the form of corn-derived (Zea mays) 



ferredoxin-NADP^ reductase (zwFNR). This enzyme 
Unks 50FD to the endogenous NADPH pool by catalyz- 
ing redox exchange between the two electron carriers 
[47]. FNR requires no maturation factors and is unaf- 
fected by O2, therefore it serves as a hydrogenase-inde- 
pendent control source of electrons for our pathway. 
Expression of zwFNR with soFD and zmSlR rescued 
growth of the cysl mutant on sulfate in aerobic selective 
media. The growth rescue required all three factors as 
well as IPTG induction of the expression plasmids. This 
result established that the ferredoxin and sulfite reduc- 
tase components of our pathway were functional and 
insulated from native metabolism under aerobic 
conditions. 

However, we found that heterologous zwFNR was not 
required to rescue growth when the selective strain was 
grown without ambient O2 (Figure 2A). While soFD and 
zwSIR expression were both still essential, this undesir- 
able background growth indicated that the native meta- 
bolic machinery could donate electrons to ferredoxin in 
anoxic conditions. We did not observe anaerobic growth 
upon expression of zwSIR alone, indicating that elec- 
trons were entering the pathway through soFD. 

Anaerobiosis effects global physiological adaptations in 
E. co/i [48]. The transition from respiratory to fermenta- 
tive growth is accompanied by a drop in cytosolic redox 
potential concomitant with an excess of reducing 
equivalents generated by glycolysis. A variety of meta- 
bolic pathways are expressed specifically in anaerobic 
conditions to dispose of electrons through electron 
acceptors including acetate, fumarate, nitrate and H2. 
Sulfite in our system represents a high potential electron 
sink, consistent with the increased tendency for elec- 
trons to enter the pathway under these conditions. 

To produce the strictest possible connection between 
the heterologous pathway and strain fitness under selec- 
tion, we sought to identify and eliminate any endogen- 
ous anaerobic electron sources. We selected 6 candidate 
interacting genes as potential nonspecific electron 
donors. Candidates were identified based on their 
homology to known ferredoxin-reducing proteins, with 
preference given to genes known to be induced 
anaerobically. 

The following six candidate genes were deleted in the 
BL21(DE3) parent strain: fpr, ydbK, her, yeaX, hcaD, 
frdB. The genes were deleted serially, in a single host 
strain, in order to expose and eliminate potentially 
redundant or epistatic interactions of the candidate 
genes with our pathway (Figure 2B). The fpr gene 
encodes an anaerobic flavodoxin-NADP* reductase [49]. 
Overexpression of ydbK, a putative pyruvate:flavodoxin 
oxidoreductase, has been shown to drive hydrogenase 
activity through ferredoxin in E. co/j [50]. The her locus 
encodes an anaerobically expressed NADH 



Barstow et al. Journal of Biological Engineering 201 1, 5:7 
http://www.jbioleng.Org/content/5/l/7 



Page 7 of 1 5 



oxidoreductase that catalyzes the reduction of the hybrid 
cluster protein Hep, an iron-sulfur protein with some 
homology to ferredoxin[51]. YeaX is a predicted oxidor- 
eductase bearing Fe-S clusters that may associate with 
the ferredoxin-like YeaW. HcaD encodes a Ferredoxin: 
NAD* reductase involved in the degradation of 3-phe- 
nylpropionate. FrdB is an Fe-S protein involved in the 
anaerobic reduction of fumarate as a terminal electron 
acceptor. Within the 4-subunit menaquinol-fumarate 
oxidoreductase complex, FrdB shuttles electrons from 
the quinone pool to the catalytic flavoprotein FrdA[52]. 

We further deleted catalytic subunits of the three 
characterized endogenous hydrogenases, hycE, hyaB and 
hybC[53], and the putative but normally silent hydroge- 
nase hyfG[54:]. All native E. coli hydrogenases are of the 
NiFe class, unrelated to the FeFe class and therefore 
unlikely to interact directly with our pathway. However, 
eliminating all native hydrogenases ensures that any H2 
production or consumption in our system could be 
attributed to the exogenous hydrogenase. 

The knockout strains were transformed with zmSIR 
and 50FD. Strains also received zwFNR as a synthetic 
electron source or an empty vector control. As 
described in the methods, growth was assayed anaerobi- 
cally overnight in selective liquid media. The results of 
these experiments are shown in Figure 2. 

The sequential knockout of candidate ferredoxin- 
interacting genes from the selection host improved the 
insulation of the test pathway from the endogenous 
redox pool. The first deletion, ^r, produced no measur- 
able effect on strain growth. The largest contribution to 
the elimination of background growth came from the 
knockout of ydbK, with further deletions only modestly 
decreasing the background growth. These results are 
consistent with the effects of individual genes observed 
by Agapakis and colleagues[55]. The mutations had no 
individual or cumulative deleterious effect on strain 
growth when zwFNR was expressed. We designate as 
EWll the final BL21(DE3)-derived strain, which bears 
the following complete genotype: E. coli B F' dcm ompT 
hsdSir^' mg') gal X{DE3) cysl fpr ydbK her yeaX hcaD 
frdB hycE hyaB hybC hyfG. All subsequent experiments 
were performed in EWll cells. 

Sulfide production by our pathway was confirmed 
spectophotometrically by the formation of methylene 
blue, as described in the methods. Wild-type BL21 cells, 
cultured anaerobically in defined medium, produce 
small amounts of sulfide: 10 |iM (± 10 at 95% confi- 
dence). No sulfide was detected in cultures of EWll 
host cells, consistent with the deletion of cysl. Expres- 
sion of soFD and zwSIR in EWll resulted in sulfide 
accumulation only to a mean level of 2 \^M. (± 1). When 
zmFNR was also expressed as a source of electrons, sul- 
fide levels increased dramatically to 200 [iM {± 30). 



Similarly, when the cflHydA hydrogenase and matura- 
tion factors were expressed, sulfide levels reached 100 
I^M (± 27). Supplying the hydrogenase with atmospheric 
H2 further raised sulfide production to 160 |iM (± 17). 
These results are consistent with the design of our path- 
way as a synthetic source of essential reduced sulfur. 

Biochemical hydrogenase O2 tolerance in situ 

We sought to initiate our selection with a wild-type 
hydrogenase with the highest possible native activity and 
O2 tolerance. We reasoned this would improve the 
probability of evolving an enzyme with properties 
exceeding those described in nature. This also would 
allow us to perform our selection in the presence of 
some O2, reducing the observed anaerobic background 
growth. Biochemical techniques allow the in vitro deter- 
mination of purified hydrogenase activity and O2 toler- 
ance[20]. But because genetic selection can be 
performed only in vivo, we assayed hydrogenase O2 tol- 
erance in cell lysates that approximate the cytosolic 
context. 

E. coli expressing a hydrogenase derived from either 
Clostridium acetobutylicum (caHydA), Clostridium sac- 
charobutylicum (csHydA), or Chlamydomonas reinhard- 
tii (crHydA), were grown to saturation in liquid culture 
under strict anaerobiosis. In each case the hydrogenase 
was coexpressed with the requisite maturation factors 
HydEF and HydG from C. reinhardtii, which are known 
to mature clostridial FeFe-hydrogenases[56]. Culture 
lysates were exposed to O2 for fixed periods of time and 
remaining hydrogenase activity was measured biochemi- 
cally, as described in the methods. 

The three hydrogenases were found to differ in both 
O2 tolerance and maximal activity levels, as shown in 
Figure 3. The anaerobic activities of the clostridial 
hydrogenases, caHydA and csHydA, were comparable to 
each other and both substantially higher than the activ- 
ity of the Chlamydomonas enzyme, crHydA. While our 
assay controlled for cell density, the in situ context of 
our system did not account for possible differences in 
expression level, maturation or folding efficiency. Repla- 
cing the C. reinhardtii maturation factors with those 
derived from C acetobutylicum yielded the same in situ 
activity for all three hydrogenases (not shown). While 
the in situ assay does not reflect biochemical specific 
activities, it measures the effective activity in E. coli 
expressing each hydrogenase, the relevant parameter for 
our genetic selection. 

We also observed substantial variation in the natural 
O2 tolerance of the three enzymes. In each case the 
inactivation by O2 could be well-described by first-order 
reaction kinetics, resulting in an exponential decrease of 
activity with time. Exposure to 0.1 atm O2 partial pres- 
sure degraded activity of caHydA with a characteristic 
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Figure 3 O2 tolerance for three hydrogenases in situ Ce Is 

expressing a hydrogenase from either C. acetobutylicum (caHydA), C. 
saccharobutylicum (csHydA) or C reinhardtii (c/HydA) were exposed 
to 10% O2 at 1 atm total pressure for the indicated times. 
Remaining activity was assessed with methyl viologen, as described 
in the methods. Three biological replicates are plotted for each time 
point. Significant differences in both maximal activity and O2 
tolerance were evident. Hydrogenase inactivation by O2 was well 
described by first order kinetics, and the best-fit exponential decay 
curves are shown. The coHydA enzyme exhibited a characteristic 
half-life of 8 (± 0.8) minutes, the csHydA enzyme 2.7 (± 0.2) minutes 
and crHydA 1.0 (± 0.3) minutes, including 95% confidence intervals. 



half-life of 8 (± 0.8) minutes, including a 95% confidence 
interval. The csHydA enzyme showed a half-life of 2.7 
(± 0.2) minutes in O2, and crHydA activity degraded 
still more rapidly, with a half-life of 1.0 (± 0.3) minutes. 
The half-life measurements are intensive biochemical 
properties of the enzymes in situ, independent of possi- 
ble differences in hydrogenase expression or maturation 
levels. 

The differences in O2 tolerance are striking. The clos- 
tridial enzymes caHydA and csHydA share 81% amino 
acid sequence similarity and nearly identical domain 
architecture. C. acetobutylicum and C. saccharobutyli- 
cum also inhabit similar strictly anaerobic ecological 
niches[57]. The algal crHydA is more divergent, only 
53% similar to caHydA, yet shares the conserved cataly- 
tic domain. Chlamydomonas, a eukaryote, exhibits a 
generally aerobic metabolism. All three enzymes in this 
experiment receive identically assembled FeFe cluster 
active sites from shared maturation factors. Yet the half- 
life of caHydA in O2 is twice that of csHydA and 8 
times that of crHydA. Because caHydA and csHydA 
showed higher activity and superior O2 tolerance, we 
chose to focus on those enzymes for further study. 

The ferredoxin-hydrogenase interaction 

The synthetic pathway we propose operates in two 
redox steps, hydrogenase to ferredoxin and ferredoxin 
to sulfite reductase. Each step must be both efficient 
and well-insulated for the overall design to be effective. 



We therefore devised independent in situ assays for 
each step of the pathway, as depicted in Figure 4A. By 
evaluating the performance of various ferredoxins in 
these assays, we sought to identify the ferredoxin best 
suited for a genetic selection. An optimal ferredoxin 
would demonstrate a robust interaction with both 
hydrogenase and sulfite reductase, while remaining insu- 
lated from nonspecific interactions with the endogenous 
redox pool. 

We first tested the hydrogenase-ferredoxin connec- 
tions, employing an independent source of electrons for 
ferredoxin. In this assay, adapted from Agapakis et al. 
[55], each ferredoxin is reduced by exogenous pyruvate- 
ferredoxin oxidoreductase derived from Desulfovibrio 
africanus (daVYOR), before transferring electrons to the 
hydrogenase. The overall activity of this pathway is mea- 
sured as H2 production in vivo, described in the 
methods. 

The direction of the hydrogenase reaction is the 
reverse of that sought for the final selection pathway, i.e. 
it is producing rather than consuming H2. FeFe- 
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Figure 4 Stepwise optimization of pathway design. A) The 

proposed pathway operates in two steps, with electrons flowing 
first from HydA to FD, then from FD to SIR. Substituting PFOR for 
SIR allows an independent test of the HydA-FD interaction. 
Substituting FNR for HydA allows an independent test of the FD-SIR 
interaction. B) PFOR-driven H2 production in vivo with a panel of 
ferredoxins. H2 production in the presence of PFOR confirms the 
ability of a ferredoxin to couple with hydrogenase. H2 production in 
the absence of PFOR suggests that a ferredoxin is reduced by a 
nonspecific endogenous source, which would tend to confound our 
genetic selection. C) FNR-rescued growth with a panel of 
ferredoxins under anaerobic conditions. Growth with FNR confirms 
the ability of a ferredoxin to couple with SIR. FNR-independent 
growth suggests that a ferredoxin is not well insulated from 
endogenous redox sources. All experiments were conducted under 
anaerobic conditions with the FW1 1 strain, as described in the 
methods. Error bars are 95% confidence intervals. 
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hydrogenases are known to operate reversibly, though 
they generally function to produce H2 in nature[58]. 
The production of H2 in this case is energetically driven 
by the highly favorable decarboxylation of pyruvate to 
acetyl-CoA[59]: 

pyruvate + CoA + 2 FDox -> acetyl-CoA + CO2 + 2 FDred + 2H* AC" - -24.8 kl mol"' 

2 FDred + 2 H* 2 FDox + H2 AC" = 3.7 kj mol"' 

Whereas in the proposed selection pathway, the con- 
sumption of H2 is instead driven by the reduction of 
highly electronegative sulfite: 

2FD0X + H2 -* 2 FDred + 2H* AC" = -3.7 kJ mol"' 

1/3 SO3 + 2 FDred + 2H' ^ 1/3 HjS -^ 2 FDox + H2O AG" = -202 kl mol"' 

H2 production from caHydA, driven by daPYOR and 
mediated by a panel of candidate ferredoxins is shown 
in Figure 4B. These data reflect both the ferredoxin- 
cflHydA interaction and the ferredoxin-i^flPFOR interac- 
tion. Because this assay does not account for variations 
in the ferredoxin-iiaPFOR interaction, it could not be 
used to infer quantitatively the efficiency of electron 
transfer between a ferredoxin and caHydA. Each ferre- 
doxin had the capacity to facilitate electron transfer to 
cflHydA, as shown by the substantial increase in H2 pro- 
duction relative to cells without exogenous ferredoxin. 
No H2 production was detectable without expression of 
a ferredoxin. Consistent with previously reported results 
[55], we observed by far the most hydrogenase activity 
with the native redox partner, caFD. 

Although expression of daPYOR significantly increased 
H2 production from each ferredoxin, it was not necessary 
for in vivo hydrogenase activity. H2 production in the 
absence of daPPOR is a measure of nonspecific electron 
flow from endogenous metabolism to a ferredoxin. Any 
such flow would tend to weaken the link between hydro- 
genase activity and host viability in a genetic selection. 
An optimal ferredoxin partner would produce high 
hydrogenase activity with daPPOR and no activity with- 
out daPYOR. We found nonspecifically driven H2 pro- 
duction to be significantly less than, yet proportional to, 
(iflPFOR-driven production for all ferredoxins. 

The ferredoxin-sulfite reductase interaction 

Optimal function of the second step in our pathway 
requires a robust and specific redox exchange of the 
mediating ferredoxin with zmSlR. We therefore used 
zmPNR as an independent source of electrons to charac- 
terize the interaction between our panel of ferredoxins 
and the zwSIR. Growth was measured in anaerobic 
selective media for strain EWll with and without 
zmFNR expression. The interactions with both zwFNR 
and zmSlR contribute to the ability of a given ferredoxin 
to facilitate growth in this system. As above, growth in 



the absence of zwFNR expression reflects the tendency 
of a ferredoxin to receive electrons nonspecifically from 
the endogenous redox pool. 

The results of growth assays from the second test 
pathway are shown in Figure 4C. In contrast to the 
results of the H2-production assay, we found that the 
clostridial ferredoxin caPD showed the worst perfor- 
mance in the sulfite reductase assay, producing no sig- 
nificant growth. Each of the other ferredoxins tested, 
50FD, zmFD and crFD, were able to effectively rescue 
the sulfide auxotrophy. In the case of the zmPD and 
crFD, we also observed significant background growth 
in the absence of zwFNR expression, indicating nonspe- 
cific interactions of these ferredoxins with native meta- 
bolism. The soPD produced no observable background 
growth under these conditions. 

None of the ferredoxins tested performed optimally in 
both the sulfite reductase and hydrogenase interaction 
assays. In both cases the differences in performance may 
be attributed to the distinction between plant-type and 
bacterial-type ferredoxins. The soFD, zmPD and crFD 
proteins all belong to the plant-type class of ferredoxins. 
These proteins carry electrons in a characteristic Fe2S2 
active center[60]. In contrast, caFD is a bacterial-type 
ferredoxin with a Fe4S4 cluster[61]. Plant and bacterial 
ferredoxins share structural and sequence homology 
[62,63] and can functionally substitute for one another 
in some cases [64]. However, it seems likely that the rela- 
tive divergence of bacterial ca¥D precludes an interac- 
tion with either zmFNR or zmSIR, both of which 
natively pair with plant-type ferredoxins. 

We also noted an apparent tendency to higher back- 
ground activities for ferredoxins with higher redox 
potentials. Reactions with electrons from the endogen- 
ous redox pool may become more thermodynamically 
favorable as the redox potential of the ferredoxin 
increases. Among the plant-type ferredoxins the soFD, 
with a redox potential of -420 mV[47], showed the low- 
est nonspecific activity in both assays. The zmPD at 
-345 mV[47] and crFD at -390 mV[65], exhibited higher 
backgrounds. This is also consistent with the theory that 
ferredoxin interactions are governed by Fe-S cluster 
redox potentials [33,66]. 

We chose soFD as a best functional compromise for 
the design requirements of our pathway. This ferredoxin 
clearly demonstrates an ability to function in both redox 
steps, and showed the least nonspecific activity in each. 
It is also one of the biochemically best characterized fer- 
redoxins and a model for ferredoxin-photosystem inter- 
actions [67]. 

O2 sensitivity of hydrogenase-rescued E. coll 

We next characterized in situ the behavior of the com- 
plete synthetic rescue pathway incorporating 50FD, 
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zmSlR and a hydrogenase electron source. In particular, 
we sought to quantify the effect of O2 on the growth of 
the selection strain in conditions as similar as possible 
to those that would be encountered in a genetic selec- 
tion. Strain EWll expressing soFD, zwSIR, and hydro- 
genase maturation factors was transformed with either 
cflHydA, csHydA, or crHydA. Cells were plated at low 
density on selective plates under custom atmospheres 
with varying O2 pressures. Growth was quantified by 
measuring the size of colonies formed after three days, 
as described in the methods. 

The dose-response relationship of O2 with various 
selection strains is depicted in Figure 5. Negative control 
strains expressing only soFD, zmSIR, and maturation 
factors showed no measurable growth under any atmo- 
sphere (Figure 5A). Positive control strains expressing 
the 02-tolerant zwFNR as an electron source showed 
robust growth under all conditions. Growth for the posi- 
tive control tended to increase with increasing O2, con- 
sistent with the energetic advantages of aerobic 
metabolism in E. coli. 

Hydrogenase-rescued strains showed 02-dependent 
growth (Figure 5B-D). For each hydrogenase, growth 
was the best in 0% O2 and decreased monotonically 
until almost no growth was detectable in 10% O2. 
Growth supported by a hydrogenase was always less 
than that observed with zwFNR, with colonies less than 
half as large forming even under strict anaerobiosis. The 



cflHydA enzyme was the most effective hydrogenase at 
supporting growth over all conditions, followed by 
csHydA and crHydA. This ordering is consistent with 
the in situ biochemical properties we determined pre- 
viously. Interestingly, we detected only weak differences 
in the 02-dependent growth profiles of the hydroge- 
nases relative to the differences in their in vivo activity 
levels and O2 tolerance. For example, caHydA shows 
roughly two-fold higher activity and an eight-fold longer 
half-life than crHydA. Yet the strain rescued with 
crHydA produced colonies only about 30% smaller and 
reached 50% growth inhibition at the same O2 level. 

Although growth in this assay was strictly dependent 
on hydrogenase expression, we found that growth did 
not require the addition of H2 to the atmosphere. Larger 
colonies were produced in the presence of H2 with all 
hydrogenases at most O2 levels, yet significant growth 
was still observed upon replacement of H2 with N2. 
Representative images illustrating the effect of H2 and 
O2 on colony size are shown in Figure 5E. H2-indepen- 
dent growth was more 02-sensitive, reaching 50% inhi- 
bition at roughly 2% ambient O2 for all hydrogenases, 
against 5% for H2-supported growth. When supplied 
with H2, strains expressing respectively c«HydA, csHydA 
and crHydA were half-maximally inhibited by O2 levels 
of 5.2% (± 0.5), 5.7% (± 0.4), and 6.1% (± 0.7), including 
95% confidence intervals. Without H2, the same level of 
inhibition was respectively reached at O2 levels of 2.5% 
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Figure 5 Hydrogenase-rescued growth is 02-dependent and H2-dependent. The genetically insulated EW1 1 strain expressing 50FD, zmSIR, 
and hydrogenase maturation factors was transformed with zmFNR (A), empty vector (A), caHydA (B), csHydA (C), or crHydA (D). Cells were 
plated at low density under custom atmospheres at varying O2 levels with and without H2. Colony sizes were measured after 3 days by imaging 
and automated analysis, as described in the methods. Representative images illustrating the effects of O2 and H2 on colony size are displayed in 
(E), with green circles indicating computationally identified colonies. Growth of strains rescued by each hydrogenase decreased monotonically 
with increasing O2 levels, becoming nearly undetectable at 10% O2. Atmospheric H2 improved growth but was not required for growth, 
suggesting an alternative source of reducing equivalents for the hydrogenase. Indicated for each curve is a best-fit sigmoid of the form area = 
a/(1+exp(b'(02-c))). Error bars are 95% confidence intervals. 
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(± 0.1), 2.2% (± 0.2), and 2.3% (± 0.3). This growth 
could not be attributed to endogenous H2 production, 
as the EWll strain lacks all native hydrogenase activity. 

Ferredoxin choice strongly affected the growth of 
the selection host in the growth-response assay. The 
substitution of zwFD for 50FD in this pathway 
resulted in much larger colony sizes by area (data not 
shown). Yet cells expressing zwFD also showed signifi- 
cant nonspecific background growth in the absence of 
zwFNR or caHydA, consistent with results shown in 
Figure 4. We chose to pursue a genetic selection only 
in strains demonstrating strictly hydrogenase-depen- 
dent growth. 

Genetic selection isolates hydrogenases with modified 
surface properties 

Variant libraries of the caHydA and csHydA enzymes 
were produced by error-prone PCR, as described in the 
methods. These libraries were used to transform EWll 
cells expressing soFD, zmSlR and the hydrogenase 
maturation factors HydEF and HydG. Transformed cells 
were grown on selective media plates for 4 days under 
10% O2, a non-permissive atmosphere for cells rescued 
by a wild-type hydrogenase. 



Colonies that appeared on the selection plates after 
incubation were restreaked twice onto selective media. 
One restreak was regrown under a selective atmosphere 
at 37°C for 4 days, the second was grown on identical 
plates in open air (-21% O2). Colonies growing in open 
air were considered false positives, while colonies grow- 
ing only in the selective atmosphere were considered 
candidates to host 02-tolerant hydrogenases. The hydro- 
genase expression plasmid was isolated from candidate 
colonies and retransformed into naive selection hosts. 
Retransformed hosts were restreaked selectively as 
before to verify plasmid-linked, hydrogenase-dependent 
growth. 

Our genetic selection isolated twenty-three hydroge- 
nases that passed all restreak and retransformation tests, 
containing a total of 110 nonsilent point mutations. 
These variants were sequenced and tested in vitro for 
hydrogenase activity as described in the methods. Of the 
twenty-three hydrogenases tested, only two retained full 
activity. A further nine hydrogenases retained some H2- 
evolution activity. The largest group of mutants, 12 of 
the 23, produced no detectable H2. The locations of the 
mutations to the caHydA and csHydA hydrogenase 
structures are shown in Figure 6. Full sequence and 
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Figure 6 Point mutations isolated by genetic selection. A) Location of 1 10 nonsilent point mutations from 21 clones. Red squares indicate 
mutations that neutralize a positively charged residue, green xs mark mutations of all other classes, vertical lines indicate stop codons. Also 
indicated are the known hydrogenase structural domains and surface-accessible domains that were inferred as described in the methods. B) The 
observed and expected frequency of mutations by charge change class. Nonsilent mutations were grouped by the charges of both the original 
and mutant amino acid. Monte Carlo simulations, described in the methods, produce the expected frequency of each charge change class. As 
indicated by an asterisk, mutations in which a positively charged residue is replaced with a neutral are overrepresented (p-value: 1.7-10"^). Error 
bars are 95% confidence intervals. 
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activity data for the mutants are provided in additional 
file 3. 

Seven of the partially functional hydrogenases retained 
complete structures, containing only point mutations. 
However, two of the partially functional hydrogenases 
were truncated by stop codons introduced into their 
coding sequences. The csHydA variant BB22 lacked 33 
residues, removing part of the catalytic domain, yet pos- 
sessed ~ 5% of wild-type activity. The cflHydA mutant 
BB05 lacked 82 residues yet still produced detectable 
H2. To our knowledge, this is the first report that such 
a severely truncated hydrogenase, lacking the small C- 
terminal subunit, remains capable of H2 production 
activity and electron transfer. 

Five of the H2-nonproducing hydrogenase mutants 
displayed only point mutations, while seven showed 
truncated catalytic domains. Interestingly, four of these 
seven mutants completely lacked the catalytic domain, 
containing only the N-terminal ferredoxin-like domains. 
The most truncated mutant that passed the selection, 
BB12, retains only 103 of 574 residues, and contains at 
most two of the four ferredoxin-like domains normally 
present. 

A recurring motif was the modification of positive 
charges, as illustrated in Figure 6. Of the 110 nonsilent 
point mutations identified by the assay, 18 altered lysine 
residues. Of these 18, 4 were mutated to arginine: 
retaining the positive charge of the lysine residue, but 
substituting a larger side chain. An additional 5 lysine 
residues were mutated to asparagine, changing the 
charge to the residue to neutral. Residue 390 was 
mutated twice, once to isoleucine and once to aspara- 
gine, in both cases neutralizing the positive charge. 
Monte Carlo simulations, as described in the methods, 
confirmed that mutations neutralizing a positively 
charged residue were statistically enriched among our 
mutants (p-value: 1.7-10'^). Because each isolated variant 
contains multiple substitutions, we were not able to 
attribute functional changes to individual point 
mutations. 

Conclusions 

The genetic selection for O2 tolerance produced no 
hydrogenases with O2 tolerance exceeding that of the 
wild-type, and only one with a comparable activity level. 
Given the observed widespread loss of H2-production 
activity, we speculate that the mutants identified by 
genetic selection confer growth by enhancing H2-inde- 
pendent electron transfer to ferredoxin. Consistent with 
this hypothesis are the results of Figure 5, indicating 
that hydrogenase, but not H2, is essential for growth 
under selection. 

We found significant enrichment of mutations neutra- 
lizing positive surface charges of the hydrogenase. 



Electrostatic forces are known to have an important role 
in the kinetics and specificity of intermolecular interac- 
tions [68,69]. Ferredoxin proteins such as those used in 
our pathway display numerous and conserved negatively 
charged surface residues, which are thought to govern 
the specific recognition of various ferredoxin redox part- 
ners [66]. Site-directed mutagenesis of lysine residues on 
the surface of Anabaena FNR was found to block inter- 
action with its native ferredoxin[70]. Our selection path- 
way employs a non-native pairing of spinach ferredoxin 
and clostridial hydrogenase, invoking suboptimally coa- 
dapted electrostatic interactions. The mutation of sur- 
face lysines may enhance ferredoxin-hydrogenase charge 
complementarity[71]. This could allow for more efficient 
electron transfer to ferredoxin, an essential activity for 
host viability. 

Our selection also isolated a number of highly trun- 
cated hydrogenase variants which were found to retain 
some function. Variants BB22 and BB05 lacked portions 
of the small C-terminal hydrogenase subunit but were 
still competent for H2 production. Similarly truncated 
putative hydrogenases have been identified in the ter- 
mite hindgut metagenome[72]. Even more severely trun- 
cated variants showed no hydrogenase activity, but were 
nevertheless capable of rescuing growth. Variant BB09, 
for example, retained only the N-terminal ferredoxin- 
like domains, which were nevertheless sufficient to facil- 
itate electron transfer to our synthetic pathway. Future 
work to identify the minimal structural elements 
required for hydrogenase function may help to structu- 
rally integrate the enzymes with designed electron flows. 

Our pathway allows 02-tolerant electron sources such 
as zwFNR to be distinguished from 02-sensitive sources 
such as cflHydA by their effects on host fitness. There- 
fore an 02-tolerant hydrogenase, once produced, could 
in principle be isolated using our selection. That no 
such hydrogenase was found suggests that other muta- 
tions exist to alter hydrogenase redox activity indepen- 
dently of the described 02-sensitive catalytic core[26]. 
Mutations of this sort may be more common than those 
specifically altering properties of the active site. The 
evolution of O2 tolerance may also require more simul- 
taneous mutations than were sampled here, or more 
extensive structural alterations. The structural features 
that optimize O2 tolerance might also change when the 
direction of hydrogenase activity favors consumption 
versus production. Future efforts may benefit from com- 
bining genetic selection with high-throughput techni- 
ques to biochemically characterize hydrogenases, 
currently in development[73]. 

The data presented in Figure 5 suggest that each 
hydrogenase can be reduced by an unknown endogen- 
ous electron source other than H2. This electron source 
is eliminated by O2, but apparently through a different 
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mechanism than that which directly inactivates the cata- 
lytic H-cluster[26]. Such a model would explain the 
apparent discrepancy between our in vivo and in vitro 
O2 tolerance assays. While Figure 2 shows different 
kinetics for the O2 inactivation of each hydrogenase in 
vitro, Figure 5 shows all three enzymes support compar- 
able 02-tolerant growth. Hydrogenase activity in our 
system might also be limited by O2 sensitivity of the 
maturation factors, rather than the mature enzymes. To 
our knowledge, potential interactions of HydEFG with 
O2 have not yet been directly examined. An in vivo 
selection system would become even more valuable in 
such a case, as mutagenesis and selection could natu- 
rally be extended to the maturation factors. 

We have shown that engineering can successfully 
insulate a synthetic electron transfer pathway from the 
endogenous E. coli redox pool. Minimizing losses to the 
cell through insulation of an artificial pathway allows 
more rational control of an engineered metabolic flux. 
We have demonstrated the use of convenient in vivo 
assays to validate isolated components of a synthetic H2 
metabolism. The results of our assays revealed trade-offs 
in the choice of pathway components, allowing compro- 
mises to meet design goals. Finally, we have successfully 
tied an essential part of cellular metabolism, the synth- 
esis of cysteine, to hydrogenase activity. By eliminating 
or disabling this activity with O2, we can halt cellular 
metabolism. These results demonstrate the utility of this 
pathway in a genetic selection for 02-tolerant hydroge- 
nases. We anticipate that future work to characterize 
the H2-independent hydrogenase activity, and to opti- 
mize ferredoxin-hydrogenase electron transfer, will allow 
for more strict selection of H2 catalysts with desired 
properties. 

Additional material 



Additional file 1: Vector sequences. Contains complete sequences for 
all vectors used in this study and described in table 1. 

Additional file 2: Media recipes. Contains exact formulations of the 
selective and induction media used in these experiments. 

Additional file 3: Mutant summaries. Contains mutation, truncation 
and activity data for all mutants isolated in this work. 
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